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^ ' Abstract. Quasar absorbers provide a powerful observational tool with 

which to probe both galaxies and the intergalactic medium up to high 

\Q ' redshift. We present a study of the evolution of the column density distri- 

bution, f{N, z), and total neutral hydrogen mass in high-column density 
quasar absorbers using data from a recent high-redshift survey for damped 
Lyman-a (DLA) and Lyman limit system (LLS) absorbers. Whilst in the 

fvq \ redshift range 2 to 3.5, ~90% of the neutral HI mass is in DLAs, we find 

that at z>3.5 this fraction drops to only 55% and that the remaining 

C ' 'missing' mass fraction of the neutral gas lies in sub-DLAs with N(HI) 

O : 1019 - 2 X 1020 cm-2. 
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1. Introduction 

Intervening absorption systems in the spectra of quasars provide a unique way to 
study early epochs and galaxy progenitors. In particular, they are not affected 
by the "redshift desert" from 1.3 ^ z ^ 2.5 where spectral emission features in 
normal galaxies do not fall in optical passbands, yet where substantial galaxy 
formation is taking place. In addition, the absorbers are selected strictly by 
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H , gas cross-section, regardless of luminosity, star formation rate, or morphology. 

Quasar absorbers are divided according to their neutral hydrogen column den- 
sities: DLAs have N(HI) > 2 x 10^'^ cm^^^ Lyman-limit systems have N(HI) 
> 1.6 X 10^^ cm^2 aj^j any system below this threshold is known as the Lyman- 
a forest. The analysis presented here is based on a sample of quasar absorbers 
found in the spectra of 66 z ~ 4 quasars (Peroux et al. 2001b) combined with 
data from the literature (Storrie-Lombardi & Wolfe 2000). 

2. Quasar Absorbers Number Density and Column Density Distri- 
bution 

The absorption lines evolution is usually described with a power law of the form: 
n{z)dz = no(l -|- z^dz, where n{z) is the observed number density of absorbers. 
For convenience, the dz is dropped and the differential number density per unit 
redshift is expressed as follow: n{z) = no(l -|- z)"^ . The observed number density 
of absorbers is the product of the space density and physical cross-section of 
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Figure 1. Number density of quasar absorbers with (from top to 
bottom) log N{HI) > 17.2, > 19.0, > 20.3 and > 21.0 07X1^"^. The 
column density range 10^^ < N{HI) < 10^^'^ cm~'^ are not direct 
observations but re-computed from the F-fit to the column density 
distribution including the expected number of LLS. The dashed line 
is a power-law fit to the number density of LLS with Zabs > 2.4. No 
absorber with log N{HI) > 21.0 cm~'^ are observed at z > 4 and the 
arrow indicates the 50% confidence upper limit. 



the absorbers which are a function of the geometry of the Universe. For no 
evolution of the properties of the individual aborbers in a A = Universe, this 
yields 7 = 1 for Qm = and 7 = 0.5 for Om = 1. We found 7 = 2.45+[J:^^ 
indicating evolution at the ~ 2a level independently of J^a/ • Figure 1 shows the 
number density of various classes of quasar absorbers. 

The column density distribution is determined as follow: 



f{N,z)dNdX 
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where n is the number of quasar absorbers observed in a column density 
bin [N, N + dN] obtained from the observation of m quasar spectra with total 
absorption distance coverage Yl^i^^i- The distance interval, dX, is used 
to correct to co-moving coordinates and thus depends on the geometry of the 
Universe. In a non-zero A-Universe: 
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The column density distribution was usually fitted by a power law over 
a large column density range N{HI) lO^'^ — 10^^ cm~^. This suggested that 
all classes of absorbers arise from the same cloud population (Tytler 1987). 
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Figure 2. Column density distribution, f{N,z), at Zabs > 3.5. The 
low column density data are Keck-HIRES observations of the Lyman- 
a forest (BR 1033-0327 and QOOOO-26, Williger et al. 1994 and Lu 
et al. 1996, respectively). The light grey bins (in the range 17.2 < log 
N{HI) < 20.3) are deduced from the fit to the observed cumulative 
number of quasar absorbers. The turn-over at the low column density 
end is incompleteness due to a combination of spectral resolution and 
signal-to-noise. The dashed and dotted lines are the two < z >~ 2.8 
power law fits from Petitjean et al. (1993) corrected for the absorber 
number density evolution with redshift and to Qa = 0.7, Om = 0.3 
cosmology. 



Assuming randomly distributed spherical isothermal halos, f{N, z) was well ap- 
proximated with a power law of slope —5/3 (Rees 1988). Nevertheless, as the 
quality of the data increased, deviations from a power law have been observed. 
In particular Petitjean et al. (1993) observed a change in slope: /5 = —1.83 for 
N{HI) < 10^^ and /3 = -1.32 above that threshold (Figure 2). 

However, LLS line profiles cannot be used to directly measure their column 
densities in the range 10^^'^ to 10^^" cm~^ because the curve of growth is de- 
generate in that interval. In our analysis, we use the expected number of LLS to 
provide a further constraint on the cumulative number of quasar absorbers and 
as clear evidence that a simple power law does not fit the observations: 



LLS, 
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where Zmin and Zmax is the redshift path along which quasar absorbers 
were searched for. We choose to fit the data with a F-distribution (a power law 
with an exponential turn-over) which was introduced by Pei & Fall (1995) and 
Storrie-Lombardi, Irwin & McMahon (1996b): 
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Figure 3. Column density distribution of quasar absorbers for var- 
ious redshift ranges. The redshift ranges are chosen to match the 
bins in the Qni+Heii plot (see Figure 4). The light grey bins (in 
the range 17.2 < log N{HI) < 20.3) are deduced from the fit to the 
observed cumulative number of quasar absorbers. The solid line is the 
F-distribution fit for z > 3.5 and the dashed lines are the fits to the 
unbinned data in the relevant redshift range. 



f{N,z) = (/,/7V,)(iV/7V,)-^e-^/^* 



(4) 



where A^ is the column density, N^, a characteristic column density and /* a 
normalising constant. Figure 3 shows the differential column density distribution 
of quasar absorbers with the F-distribution fit for various redshift ranges. The 
redshift evolution indicates that there are less high column density systems at 
high-redshift than at low redshift, confirming the earlier results from Storrie- 
Lombardi, McMahon & Irwin (1996a) and Storrie-Lombardi & Wolfe (2000). 
This suggests that we are observing the epoch of formation of DLAs. 



3. Cosmological Evolution of Neutral Gas Mass 

The mass density of absorbers can be expressed in units of the current critical 
mass density, pcrit, as: 
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Figure 4. The black circles show the neutral gas in Damped Lyman- 
a galaxies in a Oa = 0.7, Qm = 0.3 and h = 0.65 Universe (left 
panel). Vertical error bars correspond to l-o" uncertainties and the 
horizontal error bars indicate bin sizes. The light grey stars are the 
total HI+Hell including a correction for the neutral gas contained in 
sub-DLAs. The circles at low redshift are the measurements from Rao 
& Turnshek (2000), who used a method involving the observations 
of quasar spectra with known Mgll systems. The triangle at z = 
is the local HI mass measured by Natarajan & Pettini (1997). The 
squares, Qfhp, ^go and ^cetai. (Fukugita, Hogan & Peebles 1998, 
Gnedin &; Ostriker 1992 and Cole et al. 2000 respectively) are ^baryons 
in local galaxies. The right panel is for a ^m = 1 cosmology. This plot 
shows that the geometry of the Universe affects the absolute value of 

^HI+Hell with respect to the local ^baryons- 



where //, the mean molecular weight is 1.3 and mn is the hydrogen mass. 
The total HI may be estimated as: 
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If a power law is used to fit f{N,z), up to 90% of the neutral gas is in 
DLAs (Lanzetta, Wolfe & Turnshek 1995), although an artificial cut-off needs 
to be introduced at the high column density end because of the divergence of the 
integral. If instead a T-distribution is fitted to f{N, z) this removes the need to 
artificially truncate the high end column distribution and can be used to probe 
in more detail the neutral gas fraction as a function of column density and how 
this changes with redshift. The resulting VtHi+Heii is shown in Figure 4. 



4. Results and Discussion 

We find that at z>3.5 the fraction of mass in DLAs is only 55% and that the 
remaining fraction of the neutral gas mass lies in systems below this limit, in the 
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so-called "sub-DLAs" with column density 10^^ < N(HI) < 2 x lO^o cm~2 (Per- 
oux et al. 2001a). Our observations in the redshift range 2 to 5 are consistent 
with no evolution in the total amount of neutral gas. Under simple assump- 
tions of closed box evolution, his could be interpreted as indicating there is little 
gas consumption due to star formation in DLA systems in this redshift range. 
Similarly, at z > 2, Prochaska, Gawiser & Wolfe (2001) conclude that there is 
no evolution in the metallicity of DLA systems from column density-weighted 
Fe abundance measurements in DLAs (see also Savaglio 2000). At low-redshift, 

recent measurements of Qni+Heii by Rao & Turnshek (2000) at z ~ 1.65 and 
Churchill (2001) at z ~ 0.05 are difficult to reconcile with 21 cm emission obser- 
vations at z = 0. Nevertheless, the cosmological evolution of the totaZ neutral gas 
mass proves to be a powerful way of tracing galaxy formation with redshift: it 
probes the epoch of assembly of high column density systems from lower column 
density units. 

Ackno'wledgments. CP would like to thank the organising committee for 
putting together a very enjoyable meeting. 

References 



Churchill, C, 2001, ApJ, (|astro-ph/010504^ ) 



Cole, S. & the 2dFRGS team, 2000, MNRAS, ( |astro-ph/0012429| ) 
Fukugita, M., Hogan, C. & Peebles, P., 1998, ApJ, 503, 518 
Gnedin, N. & Ostriker, J., 1992, ApJ, 400, 1 
Lanzetta, K., Wolfe, A. & Turnshek, D., 1995, ApJ, 440, 435 
Lu, L., Sargent, W., Womblem D. & Takada-Hidai, M., 1996, ApJ, 472, 509 
Natarajan, P. & Pettini, M., 1997, MNRAS, 291, 28 
Pei, Y. & Fan, M., 1995, ApJ, 454, 69 

Peroux, C, Mc Mahon, R. G., Stor rie-Lombardi, L. J., & Irwin. M. J. 2001b, 
MNRAS, (|astro-ph/0107045D 



Peroux, C, Storrie-Lombardi, L. J., McMahon, R. G., Irwin. M. J. & Hook, I. 
M. 2001a, AJ, 121, 1799 

Petitjean, P., Webb, J., Ranch. M., Carswell, R., & Lanzetta, K., 1993, MNRAS, 
262, 499 

Prochaska, J., Gawiser, E. & Wolfe, A., 2001, ApJ, in press, ( |astro-ph/0101029D 

Rao, S. k Turnshek, D., 2000, ApJS, 130, 1 

Rees, M. 1988, Cambridge University Press, Vol. 107 

Savaglio, S., 2000, lAU Symposium, Vol. 204 

Storrie-Lombardi, L., Irwin, M. &McMahon, R. 1996a, MNRAS, 282, 1330 

Storrie-Lombardi, L., McMahon, R. & Irwin, M., 1996b, MNRAS, 283, L79 

Storrie-Lombardi, L. & Wolfe, A., 2000, ApJ, 543, 552 

Tytler, D., 1987, ApJ, 321, 68 

Williger, G., Baldwin, J., Carswell, R., Cooke, A., Hazard, C, Irwin, M., McMa- 
hon, R., & Storrie-Lombardi, L., 1994, ApJ, 428, 574 



